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Abstract 

We present a method to enlarge the phase space of a canonical Hamiltonian System 
in order to remove coordinate singularities arising from a nontrivial topology of the 
configuration space. This "inflation" preserves the canonical structure of the system 
and generates new constants of motion that realize the constraints. As a first illustrative 
example the spherical pendulum is inflated by embedding the sphere S 2 in the three 
dimensional Euclidean space. The main application which motivated this work is the 
derivation of a canonical singularity free Hamiltonian for the general spinning top. The 
configuration space £0(3) is diffeomorphic to the real projective space MP 3 which is 
embedded in four dimensions using homogenous coordinates. The procedure can be 
generalized to SO(n). 

1 Introduction 

One of the pillars of classical mechanics is the process of deriving equations of motion via the 
Lagrangian L. It starts with some "generalized coordinates" on the configuration space Q as 
described, e.g., in |l4], ^, |J. The Euler-Lagrange equations of motion are then obtained from 
the variational principle that J C dt be extremal on solutions of the mechanical system. For 
most configuration spaces Q there are no global singularity free coordinate systems which 
could be used as generalized coordinates in the Lagrangian. In the classical examples with 
compact configuration space (e.g., spherical pendulum, spinning top, geodesic flow on the 
ellipsoid, see e.g. |l[ [8|, 0|) global coordinates with singularities are used. This is well 
suited, e.g., in order to apply the method of separation of variables to solve the equations 
of motion. For integrable systems the coordinate singularity is typically only encountered 
by orbits with special values of the constants of motion. In the nonintegrable case the 
constants of motion are absent, while the coordinate singularity is still present. Moreover, 
it is unforseeable which orbits will encounter the singularity. Especially for the numerical 
integration these coordinate systems are therefore not advisable. 
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There are two possible cures. Since Q is a manifold there exist local coordinates ev- 
erywhere, such that we obtain Lagrangians in every chart, supplemented by the transition 
maps between different charts. This approach is certainly the most general, but it is nei- 
ther very elegant nor very simple. The second approach starts with an embedding of Q, 
dimQ = n, as a submanifold of R n+fc . In order to fulfill the constraints that define the 
submanifold Q, Lagrange multipliers are introduced, and the equations of motion are ob- 
tained by standard procedures, see e.g. 0]. Kozlov [11] introduced the notion of "excessive 
coordinates", i.e. the description of a mechanical system with more coordinates than de- 
grees of freedom. The passage from a Lagrangian with multipliers to a Hamiltonian in 



excessive coordinates is described in [11]. Our approach is different in that we start with 
the Hamiltonian in (singular) generalized coordinates and employ a transformation to a 
Hamiltonian in (nonsingular) excessive coordinates, which has all the constraints as addi- 
tional constants of motion. This "inflation transformation" is described in the first part 
and illustrated for the case of the spherical pendulum. In the second part we show how 
to apply this method to the spinning top. After this work, which first appeared in my 
thesis @, was completed, I learned that Kozlov also obtained the same Hamiltonian for 
the spinning top pi]] . Nevertheless, the method at hand gives a quite different derivation 
of these equations, which also generalizes to SO(n). 



2 Inflation of Hamiltonian Systems 

Consider a Hamiltonian with n degrees of freedom 7^(q, p). We introduce n + k new 
coordinates Q by 

q = F(Q) (1) 
c q = F(Q). 

We think of the old configuration space, parametrized by n coordinates q, as extended by k 
additional coordinates c q all of which do not show up in the Hamiltonian, i.e. we introduce 
k cyclic variables. The question now is, how to introduce corresponding new canonical 
momenta P in such a way that the desired geometric constraints F(Q) are constants of 
motion of the dynamics of a new Hamiltonian depending on (Q,P). Note that we must 
allow the transformation to have singularities at the points where the old coordinates q 
have the coordinate singularities that we are going to remove. 

This is achieved by taking (|l|) as a point transformation from Q to the trivially extended 
coordinates (q, c q ). Using the generating function 

5 = F(Q)p + F(Q)c p (2) 

we recover (||) for the coordinates by construction, while the momenta are given by 




V 9 Q 

The inverse yields the desired transformation to the new momenta: 

= (djF^V- 1 ( G(Q,P) \ 
{ dQ J { G(Q,P) J' W 

With this trick (the rest are standard canonical transformations) we can show: 
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The dynamics of the Hamiltonian System /C(Q,P) = H(F(Q), G(Q, P)) is equivalent 
to the dynamics of the Hamiltonian System 7i(q, p), i.e.: 

1. The solutions (Q(i),P(i)) of the inflated system are mapped onto solutions of the 
original system (q(i),p(£)) by (F,G). 

2. F(Q) and G(Q,P) are constants of motion of the inflated system. 

3. The canonical 1-forms are equal on solutions of the systems: p<iq = PdQ. 

These statements are clear, since the transformation is (in an extended sense) a canoni- 
cal transformation. Once we have the enlarged system, we can transform it back to the old 
one via a standard canonical transformation, with the unusual result that k coordinates 
and k momenta become cyclic. This is possible because the new Hamiltonian is singular, 
i.e. the Hessian of 7i with respect to p is degenerate. Because of this we can not pass back 
to a Lagrangian in these variables. The above three statements are comprised in the trans- 
formation of the Poisson bracket. With the notation X := (Q,P)* and x := (q, c g ,p,c p )' 
we obtain, for any function /, 



{f,K} x = (V x /)*JVx/C 
t _ 

V x f t JV x H = {f,'H} g , p . 



(5) 
(6) 
(7) 



Since 7i does not depend on c q , c p , in the last row we obtain the bracket on the original 
space. 

The Spherical Pendulum 

As a toy example we consider the spherical pendulum described by the Hamiltonian 

1 



n 



( 2 + py 

2 l P * + ri?0 



— gcos$, 



(8) 



where the singular coordinates <p) are used to parametrize the configuration space S 2 . 
We take the obvious embedding of S 2 in M 3 given by x 2 + y 2 + z 2 = 1, i.e. we take 
the coordinates r = (x, y, zf of the real M 3 in which the pendulum moves, as excessive 
coordinates: 

^ arctany/x ^ 
= arccos z jr 
\ x 2 + y 2 + z 2 j 

With the notation p 2 = x 2 + y 2 we obtain for the Jacobian 

d(F,F) 
Or 





(9) 
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(10) 
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such that the momenta are given by 







p$ 


H 







—y x 
xz/p yz/p -p 
x/2r 2 y/2r 2 z/2r 2 




(11) 



Carrying out the transformation, one obtains a singularity free Hamiltonian for the spher- 
ical pendulum in the symmetric form 



K 



{{xpy - yp x ) 2 + (xp z - zp x ) 2 + (yp z - zpyf) 



z 
r 



A short computation shows that in addition to the old constant of motion p^ = l z = xp y — 
yp x , there are geometric constants of motion r 2 and pr/2r 2 representing the constraint. 
By an appropriate choice of initial conditions we can achieve r 2 = 1 and pr = 0, so that 
the geometric constraints are the sphere and the orientation of its tangent plane in M 3 
must be the case. 



as 



Neither this coordinate system nor Hamiltonian (12) are new (according to Klein and 
Sommerfeld [10] it was first used by Hermite). As stated in the introduction, our focus is on 
the transformation procedure, which allows the application to more complicated systems 
like the spinning top. 



3 The Spinning Top in Projective Space 

As usual we take "spinning top" to be a short hand for the motion of a rigid body around 
a fixed point - with or without a potential. The configuration space of the spinning top is 
SO (3). This manifold is parametrized (again with singularities) by Euler angles. To remove 
these singularities we choose the Cayley parametrization of SO(3), which immediately 
generalizes to Q = SO(m): Consider the antisymmetric m x m matrices A with dimQ = 
n = m(m — 1) /2 independent entries £i, £2> • • • ; £n- Now introduce an additional coordinate 
£0 and map the (n + l)-tuple £ := (£o> £i> ••■) £n) onto a matrix B E SO(m) via Cayley's 
map C 

C:^^B = (e l-A)- 1 (C l + A). (13) 

Obviously every line of £ is mapped to the same B, since C(£) = C(A£) for any A G K\ {0}. 
Thus the preimage C~ 1 (B) is a line through the origin, so we can endow the £ with 
the structure of the projective space RP n . The components of £ are then interpreted as 
homogenous coordinates, and C as a diffeomorphism 

C : MP n -> SO(m). (14) 

Thus we have k = 1 for the global homogeneous coordinates, where the constraint is 
obtained from choosing a representative for every projective line on S n , by restricting £ 
to unit length. Antipodal points of S n have the same image under C. If we identify these 
points we obtain a model of SO{m) as 

S n /{±1} ~ RP n ~ SO(m). (15) 
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For m = n = 3 our equations live on S 3 ~ SU(2), the universal covering space of MP 3 ~ 
SO (3). In this case the coordinates £ are usually called Euler parameters. For the rotation 
matrix we obtain 

2£ 2 + 2i\ - e 266 - 266 266 + 266 
Bt = 72 I 2 66 + 2^6 2£ 2 + 2£ 2 - £ 2 -2^6 + 266 | • (16) 

-266 + 266 266 + 266 2$ + 2£ 2 - £ 2 

The vector (6;6>6)* is an eigenvector of B with eigenvalue 1. Thus it is the fixed axis 
of rotation described by B. The angle of rotation a is given by cos(a/2) = 6/2- The 
description of rotations by the 4-vectors £ is directly related to the Pauli spin matrices 
respectively the algebra of quaternions, see e.g. Jl^, [7|, 14, We only remark that both 



choices of "spin", and — £, by construction give the same B. 

Comparing the matrix B* with an 50(3) matrix given in Euler angles in the convention 

used in Q, we can read off the inflation transformation to obtain a Hamiltonian with £ as 
coordinates for the configuration space and constraint £ 2 = 1: 

66 — 66 /-i 7 n 

66+66 1 ) 

c2 c2 P 2 i p2 

SO ~ ?1 ~ S2 + ^3 /-, o\ 
£2 ( 18 ) 

66 + 66 



<p 


= arctan 




= arccos 




= arctan 


c q 


= 



w = Wt^^"^)- ^72^(^0-6^3) (22) 




-66 + 66 

(20) 

The old momenta given by the new ones 7r are 

Pip = 7^(6^0 - 6tti + 6^2 - 6^3) (21) 

~ " 1 

i- • tr"" ^ - 2 

Pi> = 2 ^ 37r ° + ^ 27Tl + ^ l71 " 2 + ^ o7r3 ^ ( 23 ) 
c P = £*tt/6 (24) 

Inserting this into the original Hamiltonian in Euler angles (see, e.g., j7|, [|]) the coordinate 
singularities are removed and we obtain the very symmetric Hamiltonian 

n, _ (^06 - 7Tl6 + ^26 - 7T36) 2 (^06 ~ ^16 - ^26 + 7r 36) 2 

80i 802 
N6 + vr 1 6-vr 2 6-vr 3 6) 2 + 1 {( + ^ + ^ 
80 3 6 

(-266 + 266)S 2 + (26 2 + 26 2 - H 2 )S 3 }, 

where we have used a linear potential acting on the center of mass S, and the moments of 
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inertia 0. Introducing the notation 



V £3 J 



/ a -Co ±^ 3 t6 

6 =F& -Co ±6 
V 6 ±£2 =f£l -Co 



(26) 



(the upper signs belong to £ and the lower to £) and £, z for the third row of £ we can write 
the Hamiltonian most compactly as 



1 



C 2 



(27) 



The constants of motion induced by the constraint are £ 2 and £*7r, which can be fixed to 
1 and respectively, by an appropriate choice of initial conditions. These two constants of 
motion are related to the identities = £ = and ^7r + 7r£ = 0. The angular momentum 
in body coordinates is L = £7r/2, and in space fixed coordinates $,ir /2. The reason for these 

B£. Using = £ 2 1 this gives a factorization of 

and with the 



symmetric expressions is the relation £ 
B = /£ 2 . The total angular momentum is given by 4i 2 



£ 2 7T 2 



above choice of values for the geometric constraints it just gives I 2 = 7T 2 /4. 

Note that the Hamiltonian (and the constants of motion) are invariant under a change 
of "spin", i.e. under the transformation (£, it) — ► (— £, — tt), which is induced by the trans- 
formation to the antipodal point on S 3 . Our coordinates are from T*S 3 , but the dynamics 
takes place in T*SO(3). E.g. for a periodic orbit which is non contractable in SO(3) we 
will find twice its period in 5 3 . Most notably this applies to all the relative equilibria. 

If we introduce 7 = t£ z /£. 2 f° r the coordinates of the space fixed unit vector in in- 
direction we can write the equations of motion as 



s 2^ 



IT 



2 



(28) 



Finally we remark that since l z is also a constant of motion, we can introduce new co- 
ordinates (7, L) as given above, and perform a reduction to the standard Euler-Poisson- 
equations on T*S 2 . 



4 Discussion 

We have shown how to describe a Hamiltonian system in a configuration space with higher 
dimension than the number of degrees of freedom. The resulting description with excessive 
coordinates preserves the Hamiltonian structure. This description is always useful when 
the configuration space has a nontrivial topology such that there do not exist global coor- 
dinates. Our main motivation was to obtain singularity free and Hamiltonian equations of 
motion, which are needed, e.g. for the numerical calculation of actions in integrable cases 
as described in [0, H|. The use of Euler parameters, which are used as a global coordinate 



system, has a long history. Weierstrafi [13] and Klein &; Sommerfeld [jTO used them to 



obtain explicit solutions for integrable cases. Also Whittaker [14] and Goldstein M give 



an introduction to the description of rotations with Euler parameters. However they don't 
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give a Hamiltonian in these variables. Related descriptions are also given by Kirchgraber 
and Stiefel |J for the use in perturbation theory, and in || and for numerical integra- 



tion. Let us again remark that Kozlov [11| gives the same equations, but with a different 



derivation. Our approach using Cayleys map directly generalizes to SO(n). 

5 Acknowledgement 

The author would like to thank Peter H. Richter and Andreas Wittek for inspiring discus- 
sions. This work was supported by the Deutsche Forschungsgemeinschaft. 



References 

[1] V. I. Arnold. Mathematical Methods of Classical Mechanics, volume 60 of Graduate 
Texts in Mathematics. Springer, Berlin, 1978. 

[2] H. R. Dullin. Die Energieflachen des Kowalewskaja-Kreisels. Mainz Verlag, Aachen, 
1994. Dissertation. 

[3] H. R. Dullin, M. Juhnke, and P. H. Richter. Action integrals and energy surfaces of 
the Kovalevskaya top. Bifurcation and Chaos, 4(6):1535-1562, 1994. 

[4] H. R. Dullin and A. Wittek. Efficient calculation of actions. J. Phys. A, 27:7461-7474, 
1994. 

[5] H.-D. Ebinghaus et al. Zahlen, volume 1 of Grundwissen Mathematik. Springer, Berlin, 
2 edition, 1988. 

[6] Denis J. Evans and Sohail Murad. Singularity free algorithm for molecular dynamics 
simulations of rigid polyatomics. Molec. Phys., 34(2):327-331, 1977. 

[7] Herbert Goldstein. Classical Mechanics. Addison- Wesley, Reading, MA, 2 edition, 
1980. 

[8] C. G. J. Jacobi. Vorlesungen uber Dynamik. Chelsea Publ., New York, 1969. 

[9] U. Kirchgraber and E. Stiefel. Methoden der analytischen Storungstheorie und ihre 
Anwendungen. Teubner, Stuttgart, 1978. 

[10] Felix Klein and Arnold Sommerfeld. Uber die Theorie des Kreisels. Teubner, Leipzig, 
1910. 

[11] Valerij V. Kozlov. Symmetries, Topology and Resonances in Hamiltonian dynamics. 
Springer, Berlin, 1996. 

[12] Roland Sonnenschein. An improved algorithm for molecular dynamics simulation of 
rigid molecules. J. Comp. Phys., 59:347-350, 1985. 

[13] K. WeierstraB. Vorlesungen uber Anwendung der elliptischen Funktionen. In Mathe- 
matische Werke, volume 6. Georg Olms Verlagsbuchandlung, Hildesheim, 1897. 



7 



[14] E. T. Whittaker. A Treatise on the Analytical Dynamics of Particles and Rigid Bodies. 
Cambridge University Press, Cambridge, 4 edition, 1937. 



8 



